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Abstract 

We obtain the bounds on arbitrary linear combinations of operators of dimension 6 in the Stan- 
dard Model. We consider a set of 21 flavor and CP conserving operators. Each of our 21 operators 
is tightly constrained by the standard set of electroweak measurements. We perform a fit to all 
relevant precision electroweak data and include neutrino scattering experiments, atomic parity vi- 
olation, W mass, LEPl, SLD, and LEP2 data. Our results provide an efficient way of obtaining 
bounds on weakly coupled extensions of the Standard Model. 
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I. INTRODUCTION 



Despite the enormous success of the Standard Model (SM), we are certain that the SM is 
an effective theory with a cutoff that is much smaller than the Planck scale. A lot of effort is 
being devoted to constructing and studying extensions of the SM that predict new particles 
having TeV scale masses. An integral part of this effort is constraining the parameters of 
new models using experimental data. In many cases, the constraints are obtained by directly 
computing the deviations from the SM that are induced by new particles in a specific model. 

An effective field theory approach is a two step process. First, one integrates out all new 
heavy states and obtains effective interactions involving only fields of the SM. These new 
effective higher-dimensional operators are then used to compute the deviations from the SM 
and compare with the experimental data. Following this approach one needs to make contact 
with experimental data only once — by computing the effects of higher- dimensional operators 
on different experiments and obtaining bounds on the coefficients of such operators using 
the data. Once this step is completed, one can constrain any model just by calculating the 
coefficients of new effective operators. 

The effective theory approach is by no means new, and has been applied to the elec- 
troweak data many times. Lucid explanation on applying effective Lagrangians to precision 
electroweak measurements can be found, for example, in Refs. 0, 01 . Perhaps the best known 
example of the effective approach are the so-called oblique corrections The oblique 

corrections only modify the SM gau^e boson propagators. The formalism of oblique correc- 
tions has been extended in Refs. 0, 0|- While in many models new physics contributions 
are limited to two point functions of gauge bosons, it is not always the case. 

In this article we study bounds on the coefficients of effective operators that could be 
relevant for the physics of electroweak symmetry breaking. We assume that just above the 
electroweak scale the SU{2)l x U{1)y gauge symmetry is linearly realized and therefore the 
field content involves a scalar electroweak doublet. This determines the power counting and 
our effective Lagrangian 

£ = CsM + ai Oi (1) 

contains operators of dimension 6 in addition to the SM Lagrangian. Therefore each coeffi- 
cient has the dimension of inverse mass squared and can be conveniently represented as 
ai = In our analysis, we include 21 operators — two of which correspond to the 5* and T 

parameters Our choice of operators Oi is explained in the next section. In summary, we 
choose operators that conserve U{3)^ flavor symmetry of the SM, as well as conserve CP. 
We also restrict ourselves to operators that are stringently constrained by the usual set of 
electroweak precision measurements. There is a number of dimension 6 operators that are 
flavor and CP conserving but the available data is not accurate enough to place stringent 
bounds on the corresponding Aj's. The bounds on the scales Aj of the 21 individual opera- 
tors in our basis are all about 1 TeV or higher. In contrast, the bounds on flavor conserving 
four-quark contact interactions or operators involving quarks and gluons are usually much 
lower. While it is exciting that operators suppressed by relatively low scales A are allowed, 
such operators are not very useful for constraining new models. 

Because of the spectacular agreement between the SM and precision experiments the al- 
lowed deviations from the SM are small. Therefore, when computing the effects of operators 
Oi in Eq. we can restrict ourselves to computing only the interference terms between 
CsM and the operators Oi. In other words, we work to the linear order in the coefficients 
Oj. Expansion in powers of corresponds to expansion in ^ or where E is the charac- 
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teristic energy scale of a given process and v is the Higgs vacuum expectation value (vev). 
Given that E,v <^ Aj, neglecting the terms quadratic in is a good approximation. The 
SM predictions can be computed to arbitrary accuracy since they do not depend on Oj. We 
then compare the results with the experimental data and compute the distribution as a 
function of Oj. Since our results are linear in the coefficients Oj, is quadratic in a, and 
can be written as 

= xLin + («i - ai)Mij{aj - aj), (2) 

where cij are the values of Oj that minimize x^- 

The main result of our paper is the matrix A4ij and the vector in Eq. We would like 
to stress an obvious point here — the coefficients A4ij and di are constants that we obtain by 
fitting to experimental data. Given this, admittedly, large set of numbers one can compute 
the bounds on any linear combination of operators Oj in our basis. Integrating out heavy 
particles from an extension of the SM will, in general, lead to a set of operators Oi whose 
coefficients are correlated by the underlying theory. Since we are working in the linear 
approximation in terms of the coefficients a,, the effects of correlated operators are also 
linear and the analysis of bounds on a new model remains very simple. Readers interested 
in applying our work to their favorite model can use Eq. Q and skip all the details on how 
we obtain the distribution. 

An effective Lagrangian analysis has been performed in many cases in the past. For 
examples see Refs. li ITB, H i, EE EH ■ Ref . HE is perhaps the closest to our approach, but 



it only includes the operators constrained by the Z-pole measurements and considers bounds 
on individual operators. We focus on a broader set of effective operators and use the results 
from LEPl, SLD, LEP2, as well as lower energy experiments to bound the coefficients of the 
effective operators. The most important set of precision electroweak measurements come 
from LEP and SLD. Since these experiments no longer collect data, the precision of the data 
relevant for our analysis is not going to significantly improve in the near future. 

In the next section we enumerate the operators in our basis. We further explain what 
motivates our choice of the 21 operators. In Sec. IIIH we list the experiments that we use to 
obtain the bounds on the effective operators. We outline our calculations of the interference 
terms between the SM and the additional operators in Sec. I1V[ and describe the numerical 
analysis in Sec. We summarize the results in Sec. IVll Our numerical results are presented 
in Appendix^ We show how to use these results in Appendix^ by reproducing the bounds 
on the S and T parameters. Another example is provided in Appendix O wh ere we illustrate 
our procedure on Z' gauge bosons and compare our results with Ref. We performed 
the numerical analysis using Mathematica, and we make the notebooks available online. A 
few comments about our code are contained in Appendix iDl 



II. OPERATORS 

We assume that just above the electroweak symmetry breaking scale, the effective theory 
is that of the SM with one Higgs doublet. In our effective theory the SU{2)l x f/(l)y sym- 
metry is linearly realized. The assumption of one Higgs doublet is not at all relevant since 
it is only the vev of the Higgs that enters the analysis. Since the electroweak symmetry 
is linearly realized the physical Higgs and the electroweak breaking vev are assigned mass 
dimension one. The vev appearing in the Lagrangian is always raised to positive powers. In 
case of nonlinearly realized electroweak symmetry, the expansion parameter is the momen- 
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turn divided by the vev times Air. A number of authors studied electroweak precision data 
in nonhnear reahzations of electroweak symmetry, for examples see Refs. 0, 13 • 

A complete set of dimension-6 operators consistent with the SU{3) x SU{2)l x [/(l)y, 
baryon and lepton number conservation and the linearly realized electroweak symmetry has 
been presented in Ref. 0. There are 80 operators in the basis of Ref. 0] after the leading 
order equations of motion are used to obtain an independent set of operators 1^ l^- 



We are interested in constraining models of new physics pertinent to the electroweak sym- 
metry breaking. Processes that contribute to flavor or CP violation have to be suppressed 
by scales much higher than the electroweak scale. A typical suppression for four fermion 
operators that contribute to the K — K mass difference is about 10^ TeV 0]. The bounds 
are even more stringent in the lepton sector: about 10^ TeV suppression is required for the 
electric dipole moment of the electron 17, 3]; the limits on the fi ^ e'j decay also imply 



10^ TeV suppression for the contributing operators jH. Thus it is natural to assume that 
the electroweak scale and the scales associated with flavor and CP violation are well sepa- 
rated. Processes involving the third generation are an obvious exception. Experiments have 
limited statistics and new flavor-dependent physics at TeV scale is possible. It is plausible 
that the third generation actively participates in the electroweak symmetry breaking. This 
interesting case merits a separate study but is beyond the scope of this article. We therefore 
impose U{3)^ flavor symmetry on our operators. A different U{3) acts on the left-handed 
quarks and leptons as well as on the right handed quarks and leptons. Consequently, our 
operators are unchanged when written in terms of either mass or gauge eigenstates. 

Among the 80 operators of dimension 6 listed in Ref. [H, there are 28 operators that do 
not conserve CP or flavor U{3)^, or violate both. Among the remaining 52 operators there 
are 18 operators that only involve quark and/or gluon fields. The bounds on such operators 
are poor since the precision of hadron experiments is not comparable to that of the e~^e~ 
machines. Such operators are therefore not helpful in constraining models with new physics 
at the TeV scale. 

Thus, there are 34 operators that conserve flavor and CP and contain either electroweak 
gauge bosons or leptons and perhaps some quark fields. Six of these 34 operators are not 
observable in current experiments: either they contribute to dimension 4 couplings in the 
SM Lagrangian or they involve the Higgs doublet only. 

We are finally left with 28 operators. Seven operators in this set are of the form 

Ofp = t{frD''f)F,,, (3) 

where / represents a fermion and F^^, is the field strength for the hypercharge or weak gauge 
bosons. These operators can only contribute to the Z-pole measurements. At other energies 
the interference term between the SM contribution and the contribution of operators Ofp 
vanishes since one of the contributions is real and the other is imaginary. For the same 
reason the interference term at the Z pole is suppressed by the ratio of the Z width to the 
Z mass. We therefore neglect the operators Ofp in our analysis. 

We choose the following basis for the remaining 21 operators that are the focus of our 
work. Our notation is standard: q and / represent the three families of the left-handed quark 
and lepton fields, respectively. The right handed fields are labeled u, d, and e. We omit the 
family index which is always summed over due to the flavor U{3)^ symmetry. We adopt the 
notation of Ref. |1| with minor modifications and in a few cases our operators do differ from 
Ref. by a numerical factor. 
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The operators that contain only the gauge bosons and Higgs doublets are 

OwB = = \h^DX^ (4) 

where W^^ is the SU{2) field strength, B^^ the hypercharge field strength, and h represents 
the Higgs doublet. There are 11 four-fermion operators. These are 



— 






(5) 








(6) 


Ole = 




(g7^g)(e7^e). 


(7) 


Olu = 






(8) 


Oee = 


-(e7''e)(e7/,e), Oeu 


= (e'j^e)(ujf,u), Oed = (e-f^e)(d'yf,d). 


(9) 



The operators in Eqs. (jS} and Q involve only left-handed fields, in Eqs. ((Zj) and (jHl) 2 left- 
handed and 2 right-handed, while in Eq. all right-handed fields. There are 7 operators 
containing 2 fermions that alter the couplings of fermions to the gauge bosons 

Oil = i{h^D^'h){h^l) + h.c. Oil = i{h)a''D>'h)(l-i^a''l) + h.c, (10) 

01^ = i{h^D^h){q-f,q) + h.c, O*, = i{h^a'^D^h){q^,a'^q) + h.c, (11) 

Ohu = i{h^D^h){u-i^u) + h.c, Om = i{h^D^'h){d-i^d) + h.c, (12) 

OHe = t{h^D^^h){e^^e) + \i.c.. (13) 

Finally, there is an operator that modifies the triple gauge boson interactions 

Ow = e"^" W^^Wt^W'/. (14) 

Eqs. through (|Tl|l define our basis of the 21 operators. 

We denote the coefficients in the Lagrangian in Eq. using the same indices as the 
corresponding operators, so that the effective Lagrangian is 

= ^^SM + o-wB OwB + o-hOh + ■ ■ ■ + aw Ow (15) 

Note that the first two operators in our basis, Owb and Oh, are in a one-to-one relation with 
the S and T parameters The correspondence is 

awB = T—^S, at = -2-^T, (16) 

where (h) = ^ y^^^/Q,^' fine-structure constant, s and c are the sine and cosine of 

the weak mixing angle, respectively. The V parameter is related to a dimension-8 operator 
in our power counting scheme and therefore does not appear in our analysis. 



III. EXPERIMENTS 

The three most precisely measured electroweak sector observables: a, Gp, and Mz are 
taken to be the input parameters, from which the SM gauge couplings and the Higgs vev are 
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htandard JNotation 


Measurement 


Reference 


Atomic parity 


Qw{Cs) 


Weak charge in Cs 


[21J 


violation 


QwiTl) 


TUT 1 1 ■ Til 

Weak ciiarge m il 


fool 

I22| 


DIS 


2 2 

9l,9r 


1 J-J. ■ r AT T T 

f^-nucieon scattering irom iNuieV 


Tool 

I23| 






t'^-nucieon scattering irom CDHb and CHAKM 


[24, 25J 




K 


f^-nucieon scattering irom CLl^K 


[26J 




9v ^9 A 


z/-e scattering irom uHAKM ii 


[27] 


Z-pole 


■p 


iotal Z widtn 


I2U| 







e^e iiadromc cross section at Z pole 


ronl 

I20| 




Rf\J — e,fJ',T,b,c) 


Ratios of decay rates 


ronl 






Forward-backward asymmetries 


ronl 




Sm O^ffiQFB) 


Hadronic charge asymmetry 


ronl 
[20J 




AfU = e,fi,T,b,c) 


Polarized asymmetries 


ronl 
|2(J| 


Fermion pair 


f^fU = <1,IJ',T') 


iotal cross sections tor e e j j 


ronl 
I2U| 


production at 


^Fsif = ^''^) 


Forward-backward asymmetries for e'^e'^ ff 


[20] 


LEP2 


da e/d COS 9 


Differential cross section for e^e~ — > e^e^ 


[28] 


W pair 


daw 1 d cos 9 


Differential cross section for e^e~ — > W~^W^ 


[29] 




Mw 


W mass 


[20, 30] 



TABLE I: Relevant measurements 



inferred. Predictions for experiments are computed in terms of the inputs and the coefficients 
of the new operators. The experimental quantities we use to constrain the coefficients of op- 
erators are hsted in Table HI Detailed descriptions and references for individual experiments 
can be found in many reviews, for example in Refs. and |2ol ]. 

The list of experiments in Table|l]does not include the anomalous magnetic moment of the 
muon [s^l; one of the most precisely measured electroweak quantities. The operators that 
contribute directly to {g — 2) involve left and right-handed fields and are not U{3)^ invariant. 
There are also loop contributions from operators like Owb, Ow, and many four- fermion 
operators. Such loop contributions are divergent and require introducing counterterms in 
the form of operators excluded from our analysis due to their lack of U (3)^ invariance. An 
operator analysis of contributions to the muon (g — 2) can be found in Ref. |3j]. 

For a given observable X, our prediction can be written as: 

Xth = XsM + X! '^i-^i^ (l*^) 
i 

where Xth is the prediction in the presence of additional operators, Xsm is the standard 
model prediction and X]j caXi are corrections from our new operators. In practice, the SM 
predictions are computed to the required accuracy in perturbation theory and are well known 
for all the measurements we use. Note that the corrections Xj arise in two different ways. 
First, an operator can generate a new Feynman diagram contributing to a given physical 
process. For example, a four-fermion operator Oie enters the e^e^ yU+yU^ process as a new 
diagram, in addition to the Z and 7 exchange diagrams. We call this "direct" correction. 
Second, some operators can shift the input parameters, because they add new diagrams 
to the physical processes based on which a, Gp-, and Mz are measured. Thus, the input 
parameters determined from these observables are different from their SM values. Since 
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all the other observables are calculated from these input parameters, they will inevitably 
receive indirect corrections from the shifts. We summarize the direct and indirect effects of 
our operators in Table ITTl 

Because of their high statistics, the Z-pole data and several best measured low-energy 
observables dominate the bounds on the coefficients Oj whenever such measurements con- 
strain an operator. This is the case for the operators that shift the input parameters and 
the operators of the form Oh/, which change the couplings between the Z boson and the 
fermions. The four-fermion operators do not contribute to the Z-pole measurements at the 
linear order, that is the interference term between the SM contribution and the four-fermion 
operators vanishes at the Z pole. Therefore, to constrain four-fermion operators we have to 
include the cross sections for fermion-pair production at LEP2. We also include the differen- 
tial cross sections for the W pair production to constrain the operator O^f- There are several 
operators, in addition to Ovk, that alter the cross section for W pair production. However, 
these operators are well bounded by other measurements and the W pair production does 
not contribute significantly to the bounds on their coefficients. 



Operator(s) 


shift 


Mw 


Z-pole 


DIS 


Qw 


e+e- ^ // (LEP2) 


e+e- W+W- 


OwB 


a, Mz 




V 


V 


V 


V 


V 


Oh 


Mz 














0\i 


Gf 






V 




V 




Oh Oi, 








V 




V 




Oee 












V 




Ofq,Ol^,Olu,Old 








V 


V 


V 




Oeq , Oeu j O^d 










V 


V 






Gf 




V 


V 


V 


V 


V 


oil, O he 






V 


V 


V 


V 


V 


OhU,Ohd,Of^^j,0{g 






V 


V 


V 


V 




Ow 














V 



TABLE II: Measurements influenced by different operators. The check marks, ^/, indicate "direct" 
corrections only. When an operator contributes to one of the input parameters, the corresponding 
shift of the input parameter does affect all measurements. 



IV. CALCULATIONS 

In this section we describe the computation of the effects of dimension-6 operators. Since 
a lot of work on this topic is already available in the literature, we quote the results whenever 
available. We have independently verified all the quoted results. 

We work in the linear approximation in terms of the coefficients Oj. As we indicated 
in the previous section, there are two ways that terms linear in Oj arise. First, as a result 
of additional Feynman diagrams due to the dimension-6 operators. In this case we simply 
compute the interference terms between the new operators and the tree-level contribution 
in the SM. Second, a few of our operators redefine the input parameters inferred from the 
measurements of a, Mz, and Gf- We use the tree-level SM results and expand to the linear 
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order in the deviations induced by the coefficients Oj. One of the most transparent ways of 
deahng with the shifts of the input parameters is described in Ref. 0. 

To track down all the shifts, we find it convenient to use the following parameters in the 
SM Lagrangian: e, s, and Mz] or equivalently g, g', and Mz- They are related to the input 
parameters and the new operators as 

= el{l-2vhcawB), (18) 
Ml = Mlo{l + 2v^scawB + Y('h), (19) 



y2 



-2+2aL-4, (20) 



where 



p2 1 2 

Parameters with subscripts are the values derived in the absence of any additional opera- 
tors. Eqs. (fTHj) and (fT^ are the shifts due to the S and T parameters. 

The corrections induced by S and T are given in Ref. 0], and can be easily translated to 
our notation using Eq. ()16p. Ref. 0] does not provide formulas for LEP2, but the extension 
to LEP2 is simple. The operator Owb also contains a triple gauge boson coupling, which 
contributes to the e+e^ W^W' process. We discuss this at the end of this section. 

The operators with 2 Higgs doublets and 2 fermions, Ohf, alter the couplings between 
gauge bosons and fermions. The changes to the Z-fermion couplings, gy and g^j^, are given 
in Ref. [10]. (The v"^ in Ref. jl^ is defined as one half of our value.) These changes affect 
all measurements involving a Z-exchange diagram. We have calculated the SM tree-level 
predictions in terms of arbitrary gy and g^j^, so it is easy to obtain the corrections induced 
by Ohf- In addition, the operators 0\i and change the couplings of the W boson to the 
leptons and quarks: 

g g{l + v'^a{i) (VT - leptons coupling), (22) 
g (7(1 + f ^a|jg) (ly — quarks coupling). (23) 

These changes enter, besides the shift to Gp, the calculations of the e~^e~ W^W~ cross 
section in the v exchange channel, and the i/-nucleon scattering charged current cross sec- 
tions. 

The operator 0\i shifts the value of the input parameter inferred from Gp- All other four- 
fermion operators do not contribute to the Z-pole measurements. However, they contribute 
to the low-energy measurements and LEP2 measurements. We now enumerate the effects 
of four-fermion operators on the relevant observables: 

1. The weak charges measured in atomic parity violation experiments 

Qw{Z, N) = -2[{2Z + N)Ciu + {Z + 2N)Cid]. (24) 

Giu and Cid receive corrections 

ACi„ = -^^{-al^ + a\^ + aeu + age- aiu), (25) 
^Cid = -rp^{-ai -a\ +aed + age-aid). (26) 

4GrF 
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2. i/-nucleon scattering. The 4-fermion operators affect both the neutral current and the 
charged current cross sections. The corrections to the couphngs g^, gf, g"^, gfj are 



^9leff = ^7f^[-< + (^9l - l)ay , (27) 

A^le// = ^7f^[-< + (^9i + l)«y, (28) 

^9leff = ^7f^^-< + 2^?H)' (29) 

^9ieff = ^Tf^^-^'^^ + ^^^H)- (30) 

These corrections are "effective" in the sense that they are not corrections to the Z- 
fermion couphngs and thus the formulas above only apply to the z/-nucleon scattering. 
The corrections to the measured quantities can be easily calculated from the above 
equations, for example, gj^ = {g'lY + {gtY measured at NuTeV receive the correction 

A(^i) = 2gl/\gl,^^ + 2^^A^^,,^^. (31) 



3. v-e scattering. The corrections to the coupling g'^ and g^^ are 

1 

1 



^9v%ff = ^7f^(-«« + - (32) 
A^^e// = ^7f^(-«« + + (33) 



Again, these corrections are "effective" that is only apply to the u-e scattering process. 

4. Fermion pair production at LEP2 energies. The differential cross sections in the pres- 
ence of the contact operators are given, for example, in Ref. (33l |. 



Finally, the operators Ow and Owb alter the triple gauge boson couplings. After substi- 
tuting the vev for the Higgs doublet, the two operators yield the following couplings 



AC = iawBv'gW^W-{cA^"' - sZ^"") + 6iawW-^W^\sA''^ + cZ'i). (34) 



The tree-level e"'"e~ W^W~ differential cross section is calculated in Ref. j3J| for arbitrary 
triple gauge boson couplings. Our effective couplings, Eq. (jSD), correspond to the terms 
multiplying Ky and Xy in Eq. (2.1) in Ref. ji^. To obtain the cross section we substitute 



Ak^ = — awB, (35) 



s 
v'^ s 

Akz = — —awB, (36) 
AX^ = A\z = ^aw. (37) 



where A's denote the deviations from the SM values. 
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V. TOTAL DISTRIBUTION 



In the previous section, we described how to compute the changes of observable quantities 
induced by the operators. Given theses results we calculate the total distribution. For 
non-correlated measurements, 

X = 2 ' (38) 

X ^x 

where Xg^p is the experimental value for observable X and ax is the total error both exper- 
imental and theoretical. The experimental values for the observables are obtained from the 
references cited in Table H] The input parameters and the SM predictions, except those for 
LEP2, are obtained from Ref. jl^, where the following input parameters are used: 

mmggs = 113 GeV, nitop = 176.9 GeV, a,(Mz) = 0.1213 (39) 

in addition to Mz, a, and Gp- The uncertainty in the values of the input parameters are 
incorporated as theory errors on the SM predictions and combined later with experimental 
errors. For the DIS measurements of CDHS, CHARM and CCFR, we use the SM predictions 
in Ref. |3^, but have corrected them for the small differences of input parameters |36| . 
The sensitivities of the SM predictions for the e^e" f^f if 7^ e) cross sections at 
LEP2 and the e^e~ W^W~ cross section to small changes of the input parameters 
are negligible cornpared to experimental errors, as we have verified using ZFITTER js^l 
and RacoonWW Hj. Therefore, we use the SM predictions provided in the corresponding 



experimental references. The SM prediction for the e'^e e~^e (LEP2) differential cross 
section is calculated using the program BHWIDE js^, assuming the same input parameters 
as in Eq. 

The SM predictions agree with the experimental values well, except for a significant 
discrepancy for g1 obtained by the NuTeV collaboration [23]. We include the NuTeV result 
in our calculation, but one could easily omit this result from the calculation. We also 
note that the LEP2 results for the e"^e~ qq total cross sections are larger than the 
SM predictions fairly consistently across different energies probed by LEP2. If we tried 
to constrain the coefficients of four-fermion operators with two leptons and two quarks 
using only the e"'"e~ — >• qq total cross sections, we would get relatively weak bounds on the 
coefficients aj of such operators. Weak bounds mean that the the contributions quadratic 
in the corrections in a^'s should not be neglected, contrary to what we do. However, this 
discrepancy is not supported by other data that also constrains the same operators. A 
combined fit to all observables does not yield any coefficients large enough to invalidate 
the linear approximation. Therefore we suspect this pattern is caused by a systematic error. 

Higher order terms in perturbation theory alter our tree-level calculation of the inter- 
ference terms between the SM and the additional operators. The most important effect in 
electron-positron scattering is the initial state QED radiation. In order to assess the effect 
of the radiative corrections, we compare the bounds on the 4-fermion contact operators for 
the e'^e^ ^ e) and e+e~ qq channels with the results given in Ref. ||20j. Table 

8.13. In Ref. [i^l, the e^e" qq channel is constrained using the inclusive hadronic total 
cross sections. The e^e~ — > ^ e) channel is constrained using the total cross sections 

and asymmetries for I = fi,T and assuming equal coefficients for the contact terms with 
and T. For the purpose of comparison we use the same data sets and carry out the same 
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fits. In our final results, of course, all the data is taken into account to obtain the bounds 
on contact operators. 

The comparison was carried out for different contact interactions. (The coefficients of 
4-fermion operators in Ref. 2^1 differ from our coefficients aj by a factor of An.) Several 



radiative corrections and second order terms in are considered when obtaining bounds 
on these coefficients^. We have obtained bounds for the coefficients of operators neglecting 
radiative corrections but including second order terms in a^. Except for one case, the dif- 
ferences between our bounds and the bounds in Ref. for errors on Oj are less than 20%, 
and for centra values of are less than 0.2a. When translated to the scale A, the difference 
is less than 10%, which is satisfactory for our purposes. The exception is the bound for the 
"LR" interaction for e'^e~ —>■ qq channel, which corresponds to a/„ and aid in our notation, 
in which case the errors quoted in Ref. |2(| are much more asymmetric than our estimate of 
errors. However, the central value of the coefficients differs less than 0.2a. Since the contact 
terms with two leptons and two quarks can be constrained much more stringently by the 
low energy measurements, this discrepancy cannot significantly affect our global fit. 

For e+e~ — > (/ 7^ e) channels, we have also implemented the initial-state photonic 

correction to the order 0{a), which includes initial state soft photon exponentiation and 
hard photon emission js^l- This correction is the largest radiative correction and we have 
obtained better agreement with Ref. ji^l by including it. Except for the "LR" interaction 
in the e~^e~ qq channel mentioned previously, the discrepancies between our errors and 
central values of and the results of [2^ are within 10% and O.lcx respectively. The 
remaining discrepancy likely arises from the final state radiative correction, pair production 
correction and higher order corrections, that we have not implemented. 

We can safely neglect effects that contribute at a 10% level to our estimates of the 
coefficients a^, we have nevertheless included the first order initial state QED corrections to 
e~^e~ f~ (/ 7^ e) channels in our calculation of x^- A factor of {l + ag/n) is also used 
to account for the QCD corrections for the hadronic final states j4l[. 

Eq. ()38|) must be modified to account for the correlations between different measurements. 
There are three categories of data for which the correlations between measurements cannot 
be neglected. These are the correlations between Z-pole observables j2ol | , th e experimental 
error correlations for the hadronic total cross sections at LEP2 energies [20|; and the theo- 
retical and experimental error correlations for e"'"e~ e^e" differential cross sections 4^ . 
Including correlations, 

X\a,) = Y.ma,) - ^e%) (^');.' - ^e%), (40) 

where the error matrix is related to the standard deviations ap and the correlation matrix 
Ppg as folows 

Note that often the correlations for the theoretical, statistical and systematic errors are 
different, and one should take this into account when computing the final error matrix. 
Numerical results for the distribution are presented in Appendix 1X1 

^ The radiative corrections are not mentioned in Ref. , but are described in results of individual exper- 
iments at LEP2. See, for example, Refs. 0- 



'^Iq = (^pPpq(^q- (41) 
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VI. SUMMARY 



We have obtained bounds on the coefficients of 21 dimension 6 operators in the SM. Our 
analysis is hnear in the coefficients of these operators. Therefore, the deviations from the 
SM predictions arise as interference terms between the SM and the dimension 6 operators. 
As is often the case, integrating out heavy particles leads not to just one but to several 
operators whose coefficients are related in terms of the masses and coupling constants of 
the heavy states. A global analysis of precision electroweak measurements must take into 
account all new operators induced by integrating out heavy states and account for relations 
between the coefficients of such operators. Our analysis allows obtaining bounds not just on 
each individual operator, but on their linear combinations as well. Doing so, in the linear 
approximation, does not require complicated numerical analysis, and can be done efficiently 
using our results. Of course, if the new physics contributions are "oblique" or "universal" 
only, 0,0, one does not need the whole set of 21 operators. A subset of our operators that 
only modify SM gauge boson propagators has already been considered in Refs. 0, 0] 

Our analysis is accomplished by computing the distribution as a function of the coef- 
ficients Oj. In the linear approximation, takes the form shown in Eqs. (j21) and ()A1|) . We 
concentrated on flavor and CP conserving operators. Such operators are allowed in the SM 
when suppressed by scales of the order of a few TeV. Generic flavor and CP violating opera- 
tors must be suppressed by much higher scales. This wide separation of characteristic scales 
suggests that the electroweak symmetry breaking and the flavor and CP violating sectors 
can be analyzed independently of one another. We excluded from our analysis operators 
that are not tightly constrained by the data, for example operators involving only quarks 
and gluons. Such operators are not helpful in constraining extensions of the SM. 

The bounds on the coefficients of individual operators, by which we mean that the SM 
Lagrangian is amended by only one operator at a time, can be easily obtained from Eq. 

The 1 a bound on a coefficient is ± \fM^, where Aikk indicates a diagonal element of 
Ai and is not summed over k. The fourth roots of the diagonal elements, \l M.kk^ vary from 
1.3 TeV to 17 TeV, which is a measure of how rapidly chang function of = 

What is interesting is that the eigenvalues of M. vary over a much wider range, their fourth 
roots span from 180 GeV to 21 TeV. In particular, the fourth roots of the four smallest 
eigenvalues are 180, 250, 390 and 420 GeV. This means that four linear combinations of the 
operators are much more weakly constrained than the individual operators. The emergence 
of these "weakly bounded directions" in the operator space is an interesting byproduct of 
our analysis. Of course, one can not trust the exact bounds on the "weakly bounded" 
operators. The linear analysis is not applicable when the suppression scales are so low. One 
needs to work to the quadratic order in the coefficients to reliably determine the bounds. 
However, it is clear that the bounds on such linear combinations of operators are quite weak 
and below 1 TeV. It is interesting to find out if there are heavy particles that yield one of 
the weakly constrained combination of operators when integrated out. This possibility is 
currently being investigated. 
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APPENDIX A: THE MATRIX 

Our main result can be presented in two alternative ways 

= xlnin + («i - ai)Mij{aj - CLj) = xIm + + CLiMijaj. (Al) 

In the equation above, Xmm is the minimum of in the presence of dimension 6 operators, 
and xIa/ is the value of when all coefficients a, are zero. The matrix M. is symmetric 
and positive definite. The first equation makes apparent the values of coefficients that 
minimize x^, which we call cij. The second equation is more convenient to use if only a few 
coefficients are not equal to zero. 



ai 


awB 


a-h 




Oh 






aie 


Oiqe 




(lid 




CLi 
Vi 


4.1 10^ 
1.5 10^ 


-9.3 10^ 
-23. 


-5.0 
49. 


-5.8 
76. 


-60. 
-1.1 10^ 


-6.9 
-2.4 10^ 


-0.3 
29. 


-23. 
1.4102 


-4.1102 
-36. 


-7.8102 
-68. 


7.5 
44. 


ai 




fled 


<i 






'Aq 


O-hu 




Che 


aw 




CLi 
Vi 


-5.9102 
1.0 10^ 


-6.4 10^ 
15. 


2.3 10^ 
-6.4 10^ 


19. 

-88. 


-77. 
1.0 10^ 


14. 

1.7 10^ 


-3.0102 
71. 


98. 
63. 


4.6102 
1.8102 


-4.4102 
1.0 





TABLE III: Coefficients and Vi described in Eq. HA1|) . To obtain values of Oj one needs to 
multiply the numbers in the table times 10~^ (GeV)~2 and to obtain Vi multiply times 10^ (GeV)2. 



The two equivalent sets of coefficients Sj and Vi are presented in Table IIIII The elements 
of matrix Ai are listed in Table IIVI The dimensions of these elements are easy to read off 
Eq. ()A1|) since is dimensionless and a^'s have dimension inverse mass squared. 

The numerical values of the coefficients Sj and Aiij depend on both the experimental 
values and the SM predictions. Should any of the SM input parameters change in the 
future, this affects the best fit values di in Eq. ()A1|) . but not the matrix A4. The matrix 
Aiij only depends on the sizes of errors for different measurements. Thus, Ai would change 
if experimental precision improves in the future. The central, or best fit, values di depend 
on all quantities: central values of experiments, the errors, and the SM predictions. 

At tree-level, SM predictions depend on well-measured quantities. However the strong 
coupling constant, the top mass and the Higgs mass all enter at one- loop order. The least 
known of the three is the Higgs mass and it is interesting to know how the predictions 
change as the Higgs mass is varied. It is very easy to incorporate changes of the Higgs mass 
with respect to the assumed reference value of 113 GeV in Ref which we use for the 
SM predictions. The dominant contribution from the Higgs are corrections to gauge boson 
propagators and can be incorporated by shifting the S and T parameters 0]. Higgs mass 
different than its reference value shifts the best fit value di as follows 

where only the leading logarithm in the Higgs mass is kept ^] . This is a good approximation 
as can be verified by comparing with the exact one- loop result given in Ref. j43|. 
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9.1e4 
2.4e4 
-78. 

-3.9e4 
-lAeM 
-5.5e2 

-56. 

1.3e3 
-4.0e2 
-6.9e2 

-59. 

7.8e2 

4.2e2 
-1.7e4 

5.9e4 
-1.9e3 
-9.3e3 
-6.1e2 

1.2e3 
-2.8e4 
7.7 

awB 



7.9e3 

-51. 
-1.2e4 
-1.6e2 
-1.4e2 

-9.7 
72. 
3.8 

-6.9 

-42. 

I.le2 

-83. 
-4.1e3 

1.7e4 
-1.4e3 
-4.5e3 
-6.6e2 

4.3e2 
-4.6e3 
4.5 
ah 



5.8e2 
6.7e2 

0. 

0. 
2.8e2 

0. 

0. 

0. 
5.3e2 

0. 

0. 
1.5e2 
-43. 

0. 

0. 

0. 

0. 

-I.le2 
0. 



2.2e4 
1.5e2 
5.9e2 
3.0e2 

-1.4e2 

-l.le2 
66. 
6.1e2 

-2.1e2 
1.7e2 
9.7e3 

-3.0e4 
2.7e3 
8.7e3 
1.2e3 

-8.1e2 
9.0e3 
-4.2 



a 



•II 



2.7e3 
4.6e2 
0. 

-2.7e3 
1.2e3 

1.4e3 
0. 

-1.3e3 
-1.3e3 
-5.9e2 
-7.1e2 
-2.6e3 

-49. 
-1.2e3 
-1.4e3 

4.6e2 
0. 



2.9e3 

0. 1.3e3 
-7.4e2 0. 
-2.5e2 0. 
3.3e2 0. 

0. 2.6e2 
-9.1e2 0. 
5.5e2 0. 
8.3e2 17. 
-6.6e2 -31. 
-72. 0. 
3.5e2 0. 
-4. 0. 
-1.3e2 0. 
-1.6e2 23. 
0. 0. 



2.8e3 

-1.2e3 7.1e2 

-1.4e3 5.8e2 7.8e2 

0. 0. 0. 4.8e2 

1.4e3 -4.8e2 -7.3e2 0. 8.4e2 

1.3e3 -7.3e2-6.8e2 0. 4.7e2 

3.7e2 -3.9e2 -1.6e2 1.3e2 66. 

6.6e2 -82. -3.4e2 -32. 4.9e2 

2.6e3 -1.2e3-1.4e3 0. 1.2e3 

56. -1.4e2 -36. 0. -64. 

1.2e3 -5.1e2 -6.9e2 0. 5.7e2 

1.4e3 -6.9e2 -7.2e2 0. 6.7e2 

-4.5e2 2.5e2 2.4e2 -96. -1.7e2 

0. 0. 0. 0. 0. 

O'eu 



8.8e2 
3.8e2 
47. 
1.4e3 



5.5e4 

1.5e4 6.3e4 
-6.6e3 -8.7e3 6.0e3 



aiu 



aid 



ae 



1.8e2 -2.4e4 -3.1e4 7.7e3 2.6e4 
6.7e2 -3.7e3 -4.4e3 2.2e3 4.1e3 1.4e3 
7.3e2 3.3e3 3.6e3 4.2e2 -2.9e3 1.6e2 l.le3 
-3.0e2 -2.5e4 -3.2e4 4.5e3 1.7e4 2.3e3 -2.1e3 3.2e4 
0. 6.3 -1.7 0. 0.8 0. 0. 1.4 2.6 
Ged ai, at, af„ al„ ahu and aue aw 



TABLE IV: The elements of the matrix Ai. Since it is a symmetric matrix we do not hst the redundant elements. The matrix is equal to 
the numbers listed above times lO^^(GeV)'^. We abbreviate the powers 10" as en to save space. 



T 

0.2 




0.2S 



FIG. 1: Allowed region for S and T at 90% confidence level obtained using Eq. HB1|) . 



APPENDIX B: S AND T PARAMETERS 

In our procedure, the S and T parameters correspond to awB and at as 

Ascv'^awB rj. 

b = , I = -—ah- 

a 2a 



Setting all aj, but awB and a^, to zero in Eq. ()A1|1 . we get 

X'-XU {a„s a,.) ( ) ( "™ ) + 1.5 lOV™ - 2.3 lOV 

^ ^ \ -4.8102 5.310^ ; \^t; ^ ^ 

A simultaneous fit to S and T using the above equation gives 

5 = -0.08 ±0.10, T = -0.08 ±0.10. 

The la error quoted above is obtained by projecting the Ax^ = 1 ellipse onto the correspond- 
ing axis. The plot of the 90% confidence level contour is presented in Fig. 1. Comparing 
our results with the result in Ref. |19|], Figure 10.3, shows good agreement. The results for 
S and T in Ref. |19] do not include LEP2 measurements, which indicates that LEP2 results 
do not significantly affect the bounds on the S and T parameters. 



APPENDIX C: BOUNDS ON Z' BOSONS 



Theoretical and experimental constraints on a color singlet neutral gauge boson are dis- 
cussed in Ref. l3|. The SM gauge group is extended to include a U{l)z factor with Z' its 



corresponding gauge boson. If the Z' is heavier than the electroweak breaking scale, we can 
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integrate it out and obtain the following effective Lagrangian: 



IMz, ff, ^Mz> 1 + Off! 



-E 



l 



AM' 



-glzfZHih^Di^^D,^) +h.c. 



(CI) 



where Mz' is the mass of Z', gz is the coupling constant for U{l)z, and and Zf are the 
U{l)z charges for the Higgs doublet and fermions. These charges satsify l3| 

Zi = —3Zg, Z^ = —3Zq — Zh, Zu = Zh + Zq, Z^ = Zg — Zff- (C2) 

Assuming electromagnetic strength for the Higgs-Z' and fermions-Z' couplings one obtains 

zngz = gs, ZqQz = ±gs/3. (C3) 

The authors of Ref. ^3] considered three experimental constraints on Mz'- the bound 
on the T parameter implies a bound Mz> > 0.9 TeV; the total decay width of the Z boson, 
Tz, gives a bound Mz' > 1.2 TeV for the Zggz = —gs/3 case; and the left-right asymmetry 
of the electron, Ae, gives a bound Mz' > 1.0 TeV for the Zqgz = +gs/3 case. All the above 
bounds are given at 95% confidence level. 

The T parameter corresponds to the coefficient of the first term in Eq. ()C1|) : 



T = -^zlgl (C4) 



If we use the same value as in Ref. T = —0.02 ± 0.13, we get the same bound for 
Mz'- If we consider individual measurements Tz and A^,, and use the results described in 
Section HVl we reproduce the other two bounds. 

It is interesting to compare these bounds with a global fit to all data. Using our formula 
Eq. ()A1|) . it takes little effort to obtain a constraint that incorporates simultaneously all 
operators in the effective Lagrangian ()C1|) . Using Eqs. ()C1|) and ()C2|) we obtain all the 
non-zero coefficients af. 

ah = -2z]j(3, ahf = -znZfP, a/// = -ZfZfP, (C5) 

where j3 = gz/{4:Mzr)- (Cases such as hf = hq are understood to be a^q and so on.) 
Substituting these coefficients in Eq. ()A1|) and imposing Eq. ()C3|) . we obtain as a function 
of Mz'- We then find the bounds for Mz' at 95% confidence level: 

Mz' > 2.2TeV {zqgz = y), (C6) 
Mz' > 2.4TeV {zqgz = -^), (C7) 



which are about twice as large as the bounds given in Ref. [l^ . 
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APPENDIX D: A BRIEF INTRODUCTION TO MATHEMATICA PACKAGES. 



Numerical calculation of was done using Mathematica i^l • Our code can be obtained 



at http://pantheon.yale. edu/^zh22/e'w.html or from the authors. We provide two Math- 
ematica no tebooks: ew_chi2_calculations.nb and ew_chi2_results.nb. The second notebook 
spares readers from retyping our results by giving the distribution in the form in Eq. ()A1|1 . 
It also gives the corresponding errors and the correlation matrix for the coefficients of op- 
erators. For those who want to customize our calculation to better suit their purposes, 
the first notebook contains all the inputs, formulas, and calculations. We briefiy describe 
the structure of the program below. More details are supplied in the comments inside the 
program. 

The notebook ew-chi2-calculations.nb is coded in the following order: options, input 
parameters, measurements, theoretical predictions, and the calculation of x^- 

Three options have been implemented. First option turns on or off the initial state 
radiative corrections for fermion pair production at LEP2. (The radiative corrections for 
e~^e~ e~^e~ channel have not been incorporated so far.) Second option toggles if the 
NuTeV result is included in the calculation of x^- The last one controls the presence of 
second order terms in the coefficients of four-fermion operators at LEP2, see Sec. 

The input parameters, experimental values, and SM predictions are given next. Should 
any of these numbers change in the future, one needs to modify the program accordingly. 

Next, the deviations from the SM are calculated. All formulas discussed in Sec. II VI can be 
found there. The predictions are presented as the SM values plus corrections proportional 
to Oj. These predictions are then used to calculate the distribution. We have split the 
total to track the contributions from different measurements. 
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